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I .  INTRODUCTION 


Under  sponsorship  by  the  US  Navy  Civil  Engineering  Laboratory  (CEL)  the 
Ballistic  Research  Laboratory  (BRL)  conducted  an  experimental  investigation  of 
the  structural ^response  of  metallic  beams  of  rectangular  cross-section  to  air 
blast  loading.  The  specimens  and  a  support  fixture  intended  to  provide 
simply  supported  end  conditions  were  supplied  by  CEL.  Pre-test  predictions 
of  beam  responses  were  performed  at  BRL,  principally  for  use  in  the  design  of 
instrumentation.  These  predictions  employed  static  material  property  data 
and  empirical  air  blast  parameters  in  conjunction  with  the  REPS I L  structural 
response  computer  program.  When  the  experimental  data  becEime  available  it 
was  foimd  that  the  pre-test  calculations  overpredicted  the  response  by  as  much 
as  a  factor  of  four.  This  was  cause  for  some  chagrin  since  the  response 
methodology  had  been  previously  validated  by  comparisons  with  experimental 
results  in  numerous  applications.  In  the  conclusions  to  the  Technical  Report^ 
on  this  investigation  it  was  speculated  that  the  discrepancies  were  due  to 

(a)  significant  strain-rate  effects  in  the  specimen  material  and/or 

(b)  neglect  of  frictional  and  rotatory  inertia  effects  at  the  support 
points  which  resulted  in  resisting  moments  at  the  beam  boundaries. 

Subsequent  to  publication  of  Reference  1,  dynamic  material  constitutive 
data  were  obtained  for  the  beam  material.  These  data  were  the  results  of  a 
series  of  split  Hopkinson  bar  tests  on  both  tension  and  compression  specimens. 
An  interpretation  and  an  analytical  modeling  of  these  data  are  presented  in 
Chapter  II.  Following  this  a  detailed  formulation  which  includes  the  effects 
of  rotatory  inertia  of  the  rocker  shafts  and  friction  at  bearing  interfaces  is 
given  in  Chapter  III.  Then,  in  Chapter  IV,  this  augmented  boundary  condition 
formulation  and  the  dynamic  material  constitutive  modeling  are  employed  with¬ 
in  the  REPSIL  code  for  a  reassessment  of  the  comparison  between  predictions 
and  experimental  response. 


II.  DYNAMIC  MATERIAL  PROPERTIES 

The  beam  specimens  furnished  by  CEL  were  only  identified  as  being  fabri¬ 
cated  from  3003  aluminum  alloy.  When  the  cited  discrepancies  between  experi¬ 
mental  responses  and  predictions  became  evident,  specimens  for  both  tensile 
and  compressive  split  Hopkinson  bar  tests  were  machined  at  BRL  from  unused 
specimens.  These  dynamic  test  specimens  were  forwarded  to  the  University  of 


^C.  N.  Kingery,  N.  J.  Huffington,  Jr.,  and  J.  D.  Wortman,  "Response  of  Beams 
to  Airblast  Loading,"  Ballistic  Research  Laboratory  Technical  Report 
ARBRL-TR-02369,  September  1981.  (AD  B06063SL) 

'^J.  D.  Wortman,  "RPSLID  (A  One -Dimensional  Version  of  REPSIL),"  Ballistic 
Research  Laboratory  Memorandum  Report  ARBRL-MR-03221 ,  November  1982. 

(AD  A122336) 


Dayton  Research  Institute  (UDRI)  through  the  Air  Force  Materials  Laboratory 
for  use  in  their  facility.  A  digest  of  the  material  test  data  produced  by 
UDRI^***  follows. 


Figure  1(a)  shows  representative  uniaxial  transient  values  of  engineering 
stress  0,  engineering  strain  c,  and  engineering  strain  rate  e  for  a  tensile 
Hopkinson  bar  test.  The  stress-strain  curve  in  Figure  1(b)  was  obtained  by 
eliminating  the  time  variable  between  the  a  and  e  plots  in  Figure  1(a).  If 
one  discounts  the  "noise"  induced  by  the  UDRI  data  reduction  system  it  may 
be  inferred  that  during  most  of  the  test  there  is  a  nearly  constant  value  of 
strain  rate  e  and  that,  at  least  for  this  test,  the  material  behavior  can  be 
characterized  as  elastic,  perfectly  plastic. 

Corresponding  results  for  a  representative  dynamic  compression  test  are 
shown  in  Figure  2.  For  this  test  it  may  be  seen  that  the  strain  rate  decreases 
during  the  deformation  process,  although  the  departure  from  an  average  value 
is  not  excessive.  The  stress-strain  behavior  is  approximately  elastic, 
linear  strain  hardening. 

Another  type  of  result  is  illustrated  in  Figure  3.  This  dynamic  tensile 
test  was  performed  at  a  significantly  higher  impact  velocity  (and  strain  rate) 
than  that  for  the  result  of  Figure  1.  The  large  oscillations  in  the  stress* 
versus  time  curve  are  associated  with  wave  reflections  within  the  specimen  and 
with  non-uniform  stress  distribution  associated  with  "necking."  Consequently 
the  stress-strain  curve  of  Figure  3(b)  does  not  represent  the  dependence  of 
stress  on  strain  which  would  apply  to  a  specimen  in  a  state  of  homogeneous 
uniaxial  stress. 


The  results  from  nine  Hopkinson  bar  tests  have  been  converted  to  the 
stress-strain  variables  employed  by  the  REPSIL  code;  i.e.,  Cauchy  (true) 
stress  T  versus  Almansi  strain  E,  using 


E  =  e(l  +  Js  e) 

T  =  a(l  +  e) 

(2) 

and  the  values 

thus  obtained  are  plotted  in  Figure  4. 

It  should  be  noted  that 

the  slope  of  a 

curve  in  x  versus  E  space  is  related  to 

a  a  versus  e  curve  by 

dr  _  da  a 

dE  dc  *  1+c 

(3) 

It  should  be  noted  that  the  stress  a  is  not  measured  within  the  speeim  n ; 
rather,  it  is  calculated  from  strain  measur<^mcnts  made  on  the  surroundina 
elastic  pressure  bars. 

2.4, 

Private  cormuni cations  from  d.  J.  bless  and  A.  C%illita ,  PDRl  to  N.  J. 
Huffington,  Jr. ,  BRL. 

4 

S.  J.  Bless,  A.  Challzta,  and  A.  K.  Rajendran,  "Ifinami-c  Tensil<'  Tos1  Results 
for  Several  Metals,"  Air  For<fe  Material  Lahoratoru  Revert.  AF]rJAl,-TR-y‘'’-'U)rR , 
April  1982. 
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Figure  2.  Results  from  compression  Hopkinson  bar  test  ffHB057. 


•V 
'V  -- 


>'.s\v'.'sV*V..  . 


(iM)  ssBais 


*>dW  SS3dlS 


NlVdiS 


ss3ais 


Since  many,  but  not  all,  of  the  tensile  Hopkinson  bar  tests  had  the  form  of 
Figure  1(b)  with  da/de  =  0  and  e  <  <  1,  these  tensile  x  versus  E  curves  have  a 
linear  strain  hardening  form  with  slope  approximately  equal  to  the  flow  stress 
The  strain  rate  labels  on  Figure  4  are  average  values;  in  some  tests  the 
variation  of  E  was  considerable. 

Obviously,  there  is  a  disconcerting  amount  of  scatter  associated  with  the 
tensile  Hopkinson  bar  tests.  It  must  be  appreciated  this  test  is  more  open 
to  question  than  the  compression  Hopkinson  bar  test  since  it  is  difficult  to 
achieve  uniform  stressing  in  a  very  short  specimen  through  threaded  ends.  The 
dashed  curves  in  Figure  4  were  from  tests  which  were  regarded  as  less  reliable 

.  _1 

or  were  provided  belatedly.  The  curve  labeled  E  =  1480s  corresponds  to  the 
data  of  Figure  3,  where  necking  of  the  specimen  occurred.  The  curve  labeled 
•  -1 

E  =  369s  was  derived  from  a  test  which  exhibited  a  significant  strain 
hardening  slop<  da/de.  Only  one  compression  Hopkinson  curve  is  shown  in 
Figure  4  since  the  other  compression  tests  were  conducted  at  a  much  higher 
impact  velocity  and  provided  no  useful  data  in  the  range  of  E  shown.  These 
higher  velocity  tests  did  reveal  a  monotonic  increase  in  the  magnitude  of  the 
flow  stress  tending  toward  a  saturation  limit. 

Figure  4  also  includes  a  curve  labeled  E  =  0  which  corresponds  to  the 
static  stress-strain  data  published  in  Appendix  A  of  Reference  1.  From  an 
examination  of  this  figure  it  may  be  concluded  that,  in  spite  of  the  scatter 
in  the  dynamic  tensile  data,  the  beam  specimen  material  is  significantly  rate- 
sensitive,  entailing  an  increase  in  the  dynamic  flow  stress  of  50%  to  100% 
over  the  static  value. 

In  order  for  dynamic  material  property  data  to  be  useful  in  response 
analysis,  it  is  necessary  to  incorporate  these  data  into  a  mathematical  rate- 
dependent  constitutive  model .  In  Reference  1  a  single  post-test  response 
calculation  was  made  using  a  power  law  strain-rate  dependent  model  which  was 
available  within  the  REPSIL  code: 

"d  -  "s  *  ll 

where 

T  =  static  uniaxial  flow  stress 
s 

Tj  =  dynamic  uniaxial  flow  stress 

D,p  =  parameters  available  for  fitting  experimental  data. 

At  that  time  no  dynamic  data  were  available  for  the  beam  material  so  published 
values  of  D  =  6500s”^  and  p  =  4  for  6061-T6  aluminum  alloy  were  used  in  the 

'^P.  S.  Symonds,  "Viscoplastic  Behavior  in  Response  of  Structures  to  Dynamic 
Loading, "  in  Behavior  of  Materials  Under  Dynamic  Loading,  N.  J.  Buffington, 
Jr.,  Editor,  ASMS,  196h. 


response  calculation.  Since  6061-T6  aluminum  is  known  to  have  only  a  slight 
rate-sensitivity  it  is  not  surprising  that  this  calculated  response  differed 
by  only  a  small  amount  from  the  response  when  rate  effects  were  omitted. 


The  model  represented  by  Equation  (4)  has  a  serious  disadvantage  for 
numerical  calculations.  Because 


for  P  >  1, 


E  -  0 


Equation  (4)  requires  large  jumps  in  the  flow  stress  when  E  becomes  small  or 
changes  sign.  These  jumps  inject  non-physical  oscillations  in  time-marching 
solutions.  Inasmuch  as  the  physical  dependence  of  t  on  E  is  believed  to  be 
less  abrupt  near  E  =  0  and  to  approach  a  saturation  value  for  E  large  it  was 
decided  to  postulate  that  the  dynamic  inelastic  constitutive  function  is 

Tj(E,E)  =  t^(E)  +  K  tanh  BE  (6) 

where  K  and  6  are  parameters  which  may  be  adjusted  to  provide  correlation  with 
dynamic  material  tests. 

The  results  from  both  the  Hopkinson  bar  tests  and  the  static  tests  have 
been  replotted  as  versus  E  in  Figure  5  in  order  to  facilitate  a  graphical 

fitting  of  the  surface  defined  by  Equation  (6) .  Each  dynamic  test  covers  a 
range  of  strains  E  (and  of  strain  rate  E)  but  the  results  are  too  eratic  to 
permit  construction  of  curves  of  E  =  constant.  Values  of  K  and  B  were  deter¬ 
mined  initially  when  only  the  "open-loop*'  data  were  available  (Fit  I) .  When 
the  additional  data  shown  in  the  "hatched  loops"  was  provided  a  second  eval¬ 
uation  of  the  parameters  was  made  (Fit  Ilj.  Both  sets  of  parameters  are  list¬ 
ed  in  Table  1. 

TABLE  1.  STRAIN-RATE  PARAMETERS 


Fit 

K(MPaj 

K(psi) 

e(s) 

I 

64.8 

9392 

0.008069 

II 

126.7 

18370 

0.003604 

Curves  corresponding  to  these  fits  are  indicated  in  Figure  5  for  t^(0.08). 

For  other  values  of  these  curves  would  shift  vertically. 

The  matter  of  implementing  the  constitutive  function  of  Equation  (6)  in 
a  response  code  which  employs  the  mechanical  sublayer  model  has  been  investi¬ 
gated.  In  this  model 


the  total  stress  at  a  material  point  is  obtained  as  a  weighted  average 
of  the  stresses  in  m  hypothetical  elastic,  perfectly  plastic  sublayers,  all 
of  which  are  at  the  same  strain  E.  Each  sublayer  has  a  distinct  yield  stress 
The  are  the  appropriate  weighting  factors.  If  the  rate-dependent 
stresses  are  to  be  determined  using  a  similar  expression 


t^(E,E)  =  ^  St\(E,E) 


with  the  same  weighting  factors  Cj^  then  it  is  necessary  to  take  the  yield 
stress  of  the  first  sublayer  to  be 


^01^^^  ^01  Cj 


and  the  remaining  sublayer  yield  stresses  to  be 


^Ok  "^Ok 


k  >  1 


III.  FRICTION  AND  ROTATORY  INERTIA  MODELING 

Since  it  was  desired  that  the  beam  specimens  be  tested  with  simply 
supported  ends,  CEL  had  designed  and  fabricated  rocker  shafts  which  supported 
the  ends  of  the  test  specimens  and  which  in  turn  were  supported  in  split 
bearing  blocks.  Had  these  shafts  been  massless  and  had  all  rubbing  surfaces 
been  frictionless  this  design  should  have  provided  a  good  simulation  of  the 
desired  boundary  conditions  (mathematically  specified  as  no  transverse  dis¬ 
placement  and  no  bending  moment  at  either  end  of  the  beam) .  However,  in  view 
of  the  cited  discrepancies  between  predictions  for  response  of  simply  supported 
beams  and  the  experimental  data  it  was  decided  that  the  effects  of  imperfect 
boundary  condition  simulation  should  be  evaluated.  In  an  effort  to  avoid 
large  membrane  forces  in  the  beams  when  finite  deflections  occur,  the  test 
fixture  had  been  designed  to  permit  one  end  of  the  beam  of  slide  through  its 
rocker  shaft.  Owing  to  this  lack  of  symmetry  of  the  boundary  conditions  a 
separate  formulation  was  developed  for  each  boundary. 

A.  Formulation  for  the  Pinned  End  Condition 


One  end  of  each  beam  specimen  was  bolted  to  its  rocker  shaft  with  a  pair 
of  force  transducers  sandwiched  between  the  lower  surface  of  the  beam  and  the 
attachment  surface  of  the  shaft.  This  end  of  the  beam  could  only  rotate,  no 
translation  being  permitted  in  either  the  transverse  or  longitudinal  direction 
A  drawing  showing  details  of  the  pinned  end  rocker  shaft  is  included  as 
Figure  6.  Free  body  diagrams  showing  the  forces  and  moments  acting  on  the 
rocker  shaft  and  on  the  end  segment  of  the  beam  are  presented  in  Figure  7. 

The  symbols  have  the  following  definitions: 


(a)  Rocker  Shaft 


(b)  Beam  Segment 


Figure  7.  Free  body  diagrams  at  pinned  end. 
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6  =  angle  of  rotation  of  rocker  shaft  and  of  normal  to  beam 
reference  surface  at  x  =  0. 

p(t)  =  blast  pressure  loading 

Pq  =  membrane  force  in  beam  at  x  =  0. 

Sq  =  transverse  shear  force  in  beam  at  x  =  0. 

Mq  =  bending  moment  in  beam  at  x  =  0* 

Rq  =  reactive  force  of  bearing  blocks  on  rocker  shaft. 

=  frictional  moment  acting  at  bearing  interface. 

A  =  projected  area  of  rocker  shaft  surface  on  which  blast 
pressure  acts. 

e  -  radius  of  bearing 

y  =  coefficient  of  Coulomb  friction. 

The  equilibrium  equations  for  the  rocker  shaft  are 

Pq  cos  0  +  Sq  sin  e  +  =  0  (11] 

-Pq  sin  0  +  Sq  cos  0  +  R^^  -  pA  =  0  (12] 

vdiere  R^^,  R^^  are  the  x-  and  y-  components  of  R^.  These  components  can  be 
combined  to  obtain 


%  -  V'OX  *  “Oy 

(13] 

=  ^Pq  +  Sq  +  2pA  (Pq  sin  0  -  Sq  cos  0)  ♦  p^A^j^ 

The  frictional  moment  can  then  be  determined  using 
=  ye  Rq  C(0) 

=  ye  C(0)|^Pq  +  Sq  +  p^A^  +  2pA  (Pp  sin  0  -  Sp  cos  0)j*^ 


The  small  ’’overhang"  of  the  beam  is  neglected  except  for  its  contribution  to 
the  moment  of  inertia  of  the  rocker  assembly. 


where  C(d)  is  the  nonlinear  function  which  may  be  represented  graphically  as 
indicated  and  which  presently  will  be  formulated  analytically. 

c(6)  I 


The  rotational  equation  of  motion  from  the  free  body  diagram  of  Figure 
7(a)  is 


~  ’^fO  ^0® 


where  Iq  is  the  mass  moment  of  inertia  of  the  material  which  rotates  with  the 

angular  acceleration  6.  This  ordinary  differential  equation  may  be  replaced 

2 

by  the  customary  0 (At  )  central  difference  approximation  consistent  with  the 
time  integration  algorithm  employed  in  the  REPSIL  code.  Letting 


A0*  =  -  e(t.) 


(16) 


Ae'  =  0(t.)  -  e(t._j) 

where 


(17) 


(18) 


(19) 


By  use  of  Equation  (19)  a  discrete  numerical  solution  for  0(t)  can  be  obtained 

concurrently  with  the  marching  out  of  the  solution  for  the  beam  response.  The 

initial  conditions  for  0(t)  are  0(0)  =  0  and  0(0)  =  0.  In  utilizing  Equation 

(14)  to  obtain  the  value  of  M__  .  to  employ  in  Equation  (19),  one  must  detcr- 

,  t  u ,  j 

mine  the  value  of  C(0)  in  accordance  with  one  of  the  following  three 
possibilities: 

1.  0  >  0  or  (0  =  0,0  >  0)  for  which  C(0)  =  +1 

2.  0  =  0,  0  =  0  for  which  -1  <  C(0)  <  1  and 
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3.  6  <  0  or  (0  =  0,0  <  0)  for  which  C(0)  =  -1 


Although  somewhat  awkward  in  appearance,  these  tests  to  determine  which  regime 
is  operative  can  be  readily  programmed  into  the  time  marching  algorithm  of 
the  response  code. 

B.  Formulation  for  the  Sliding-Hinged  End  Condition 

The  support  condition  at  the  other  end  of  the  beam  is  similar  to  that  at 
the  pinned  end  in  that  the  rocker  shaft  (see  Figure  8)  does  not  translate  but 
is  complicated  by  the  provision  for  the  beam  to  slide  through  the  slot  in  the 
rocker  shaft  as  the  lateral  displacement  of  the  beam  varies.  The  formulation 
of  the  boundary  conditions  at  this  end  of  the  beam  entails  the  introduction 
of  two  additional  dependent  variables: 

jIi  =  angle  of  rotation  of  the  rocker  shaft 

u  =  displacement  of  point  in  beam  initially  on  axis  of  rocker  shaft 
parallel  to  sliding  interface 

Other  symbols  employed  in  this  formulation  are: 

S,,  =  transverse  shear  force  in  beam  at  x  =  N. 

N 

F^  =  frictional  force  at  sliding  interface 

Hj.  =  moment  transmitted  between  beam  and  rocker  shaft 

=  frictional  moment  acting  at  bearing  interface 

=  reactive  force  of  bearing  blocks  on  rocker  shaft 

B  =  projected  area  of  rocker  shaft  surface  upon  which  blast  pressure 
acts 

m  =  mass  of  beam  considered  in  free  body  diagram,  including  overhang 
material 

h  =  depth  of  beam 

Referring  to  Figure  9(a),  we  may  write  the  equation  of  motion  for  the 
direction  parallel  to  the  sliding  surface  as 


-P  +  F^  =  mu 


where  it  is  tacitly  assumed  that  the  angular  velocity  of  the  rocker  shaft  is 
never  large  enough  that  a  distinction  between  relative  and  absolute  accelera¬ 
tion  of  m  is  necessary.  Defining  incremental  changes  in  u  in  the  same  manner 
as  for  0,  Equation  (20)  may  be  expressed  in  finite  difference  form  as 
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i 


Au"  =  AU-  + 

m 


(-P  +  ViS^C(u)) 


Since  the  REPSIL  code  employs  Cartesian  components  of  displacement  increments, 
these  can  be  obtained  for  the  mesh  location  initially  on  the  rocker  shaft  axis 
by  use  of 

Au^  =  -Au  cos  4)  (22) 

AUy  =  -Au  sin  ij;  (23) 

Also  with  respect  to  Figure  9(a)  we  may  write  the  moment  equation 

M  -  Nij.  +  Se  +  I  =  Igij;  (24) 

where  I_  is  the  centroidal  moment  of  inertia  of  m. 

D 

For  the  rocker  shaft  (see  Figure  9(b))  we  obtain  the  following  equations  of 
equilibrium  and  of  motion; 

-F_  cos  4;  +  S.,  sin  4)  +  R,,  =0  (25) 

t  N  Nx 

-F^  sin  4)  -  cos  4)  +  R^^^  -  pB  =  0  (26) 

«T  -  -  I  ff  =  's  ♦  <”> 

Making  use  of  Equations  (25)  and  (26)  the  bearing  reaction  is  obtained: 

'S)  V'^Nx  ‘^Ny 

(28) 

=  [fJ  +  sJ  +  p^B^  +  2pB  (F^  sin  4-  +  cos  4^)]’"^ 

The  frictional  moment  at  the  bearing  is  then 
MfN  =  ve  C(4-) 

(29) 

=  ye  C(4i)  [F^  +  +  p  B  +  2pB  (F^  sin  ip  +  cos  4')]"" 

Adding  Equation  (24)  to  Equation  (27),  one  obtains 

M  -  +  Se  =  Jj^ip  (50) 

where  I„  =  I  +  I„  is  the  mass  moment  of  inertia  of  the  shaft  and  the  contained 
N  s  B 


Equation  (30)  may  be  expressed  as  the  difference  equation 


Alp  =  Alp  + 


(At)‘ 


(M 


N.j 


-  M 


fNJ 


e  S.,  .) 

N,)-" 


(31) 


The  formulation  of  the  boundary  conditions  for  the  sliding-hinged  end 
contains  the  two  Coulomb  functions  C^u) ,  C(i)  which  must  be  evaluated  by  tests 
analogous  to  those  prescribed  for  C(0) ;  i.e.,  the  possibilities  of  forward 
slip,  no  slip,  or  reverse  slip  must  be  considered  at  each  sliding  interface. 

The  foregoing  boundary  condition  formulations,  which  take  account  of 
boundary  friction  and  the  rotatory  inertias  of  boundary  masses,  have  been 
incorporated  into  the  REPSIL  code  as  a  temporary  modification  in  order  to 
perform  an  evaluation  of  the  contribution  of  each  of  these  effects. 


IV.  APPLICATION  OF  THE  GENERALIZED  RESPONSE  ANALYSIS 

All  calculations  to  be  reported  in  sequel  pertain  to  the  case  identified 
as  Shot  Number  5  in  Reference  1.  The  beam  specimen  had  a  rectangular  cross 
section,  5.08  cm  (2  in.)  wide  by  1.11  cm  (7/16  in.)  deep,  and  a  span  of  30.48  cm 
(12  in.);  it  was  fabricated  from  3003  aluminum  alloy.  Static  tensile  stress- 
strain  data  for  this  material  were  reported  in  Figure  A-2  of  Reference  1.  The 
blast  loading,  assumed  uniformly  distributed  spatially,  was  the  experimentally 
recorded  pressure  history  shown  in  Figure  48  of  the  same  reference.  Other 
input  parameters  used  in  the  calculations  are: 

A  =  2419  mm^  (3.75  in^) 

B  =  1452  mm^  (2.25  in^) 

C  =  15.1  mm  (0.594  in) 

Iq  =  0.15325  g-m^  (0.0013564  lb*in*s^) 

=  0.05360  g-m^  (0.0004744  lb*in*s^) 

m  =  0,0303  kg  (0.0001728  Ib-s^/in) 

The  computer  runs  of  the  REPSIL  code  which  were  made  for  correlation 
with  the  experimental  data  for  Shot  Number  5  are  listed  in  Table  2.  A  com¬ 
parison  of  transient  raid-span  deflections  for  each  run  is  shown  in  Figure  10 
as  well  as  the  experimental  final  value  of  this  quantity  (unfortunately,  the 
transient  experimental  data  were  lost).  Damping  was  initiated  in  each  run 
after  the  maximum  displacement  in  order  to  obtain  a  predicted  final  deflection, 
which  is  listed  in  Table  2.  Runs  A,  B,  and  C  were  reported  in  Reference  1  and 
are  included  for  completeness.  It  is  clear  that  use  of  simply  supported  end 
conditions  (run  A)  results  in  a  gross  overprediction  of  the  response.  On  the 
other  hand,  runs  B  and  C  which  were  made  using  fixed  end  conditions  effec¬ 
tively  bracket  the  experimental  final  deflection.  However,  run  C  employed 
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the  power  law  strain-rate  model  of  lujuation  (4)  with  parameters  for  a  different 
aluminum  alloy.  When  dynamic  data  for  300.1  aluminum  alloy  became  available 
the  parameters  of  I'it  1  given  in  Table  1  were  used  with  the  hyperbolic  tangent 
strain-rate  model  of  liquation  (6)  to  make  run  D,  also  with  fixed  boundaries. 
This  run  provided  excellent  agreement  with  the  experimental  final  deflection. 
However,  it  was  realized  that  the  physical  response  did  involve  rotation  of 
the  rocker  shafts  and  slippage  of  one  end  of  the  beam  through  the  slotted 
rocker  shaft  so  a  more  rational  response  modeling  was  sought. 

The  remaining  REPSIL  runs  employed  the  boundary  condition  formulations 
developed  in  Chapter  III.  Runs  E  and  F  were  made  using  static  material  prop¬ 
erties  (no  strain-rate  dependence  in  the  constitutive  function)  and  illustrate 
the  dominant  effect  of  the  friction  coefficient  n  on  the  response.  Run  C  was 
then  made  using  the  same  value  of  n  as  Run  F  but  employing  the  hyperbolic 
tangent  strain-rate  model  and  the  data  of  Fit  1.  As  may  be  seen  from  Figure 
10,  the  inclusion  of  strain-rate  dependence  had  a  relatively  small  effect  in 
this  instance,  amounting  to  a  further  decrease  in  the  final  deflection  of  12% 
of  the  experimental  value.  When  Fit  II  was  effected.  Run  H  was  made  using 
these  material  properties  and  the  same  value  of  n;  the  results  were  essentially 
identical  to  those  of  run  G.  This  somewhat  surprising  indifference  to  the 
strain-rate  data  fitting  has  the  following  explanation.  For  Shot  No.  .S  the 

maximum  strain  rate  calculated  was  only  120.8s  ^  at  one  location  in  the  beam 
and  for  most  of  the  response  the  strain-rate  was  much  less  throughout  the 
beam.  Referring  to  Figure  5,  one  sees  that  there  is  no  significant  difference 
between  Fits  I  and  II  for  such  small  strain-rates.  Conversely,  for  stronger 
blasts  which  induce  higher  strain-rates  one  could  expect  that  the  difference 
between  solutions  using  Fits  I  and  II  would  be  significant. 

It  is  instructive  to  examine  certain  other  results  produced  by  the  REPSII, 
program  for  Run  G  (to  which  results  for  other  runs  are  qualitatively  similar). 
Figure  11  shows  the  variation  of  the  bending  moment  and  the  frictional 

moment  at  the  pinned  end  while  Figure  12  displays  the  transient  rotation 

of  the  rocker  shaft  at  the  same  location.  It  may  be  seen  that  no  rotation 
occurs  until  about  0.23  ms  when  is  no  longer  able  to  match  Subse¬ 

quently  the  end  of  the  beam  rotates  until  about  1.81  ms  when  the  shaft  "locks 
up"  (M^-q  =  Mq)  and  0  remains  constant  until  the  reversed  bending  moment  be¬ 
comes  large  enough  to  overcome  the  frictional  resistance.  Another  brief  period 
of  ankylosis  starts  at  about  4.13  ms.  The  predicted  transient  behavior  of 
the  transverse  shear  force  at  the  pinned  end  of  the  beam  is  presented  in 
Figure  13.  It  should  be  observed  that  the  maximum  magnitude  of  this  shear 
reaction  occurs  before  the  end  of  the  beam  starts  to  rotate.  However,  prior 
to  this  the  material  at  this  location  experiences  inelastic  deformation.  \o 
explanation  for  the  discrepancy  between  predicted  values  of  this  quantity  and 
the  experiment  value  is  feasible  since  the  experimental  value  is  .30",  greater 
than  that  which  would  bo  expected  if  the  bounilary  condition  were  total  fixit\'. 

It  was  desired  to  assess  the  importance  of  including  the  rotatory  iiu-itia 
of  the  end  masses  in  the  response  formulation.  This  cannot  he  .iccomp 1 i shed 
by  simply  making  a  run  with  the  mass  moments  of  inertia  and  set  c(|ual  to 

zero  (see  l-ciuations  (19)  and  (31)).  Rather  than  program  a  special  formulati(m 


2.0 

TIME  (ms) 


for  this  case,  the  available  program  was  used  with  the  values  of  and 

set  to  a  small  fraction  of  their  actual  values.  This  was  done  in  Run  1  where 
the  fraction  employed  was  1/10.  The  results  should  be  compared  with  Run  G 
since  all  other  parameters  were  the  same  for  these  two  runs.  It  may  be  seen 
from  Figure  10  that  the  form  of  the  response  curve  is  somewhat  altered  for 
Run  I.  The  final  displacement  is  about  lb%  (of  the  experimental  value) 
greater  when  the  rotatory  inertia  effect  is  minimized. 

Since  none  of  the  cases  examined  for  a  =  1  resulted  in  a  final  displace¬ 
ment  approximating  the  experimental  value  it  was  decided  to  explore  further 
the  effect  of  varying  this  parameter.  This  was  done  in  Runs  J  and  K  using 
a  =  2  and  3  respectively.  Run  K  provided  a  very  satisfactory  correlation 
with  respect  to  final  deflections.  In  fact,  further  increase  in  u  would  have 
little  effect  since,  even  for  Run  K,  the  end  masses  were  "locked"  during  most 
of  the  run.  The  use  of  such  large  values  of  friction  coefficient  may  seem 
inappropriate  but  there  is  essentially  no  data  on  Coulomb  friction  under 
conditions  where  both  the  driving  force  and  the  normal  pressure  are  rapidly 
changing. 

It  is  observed  that  strains,  both  experimental  and  analytical,  vary  di¬ 
rectly  with  the  transverse  displacement  amplitude.  Thus,  the  best  correlations 
with  experimental  strain  data  occur  for  runs  which  best  match  the  experimental 
final  deflection. 


V.  CONCLUSIONS 

The  foregoing  analytical  formulations  for  strain-rate  effects  plus  fric¬ 
tion  and  rotatory  inertia  at  supports  have  permitted  a  more  rational  inter¬ 
pretation  of  experimental  data  from  blast  loaded  beam  tests  reported  in 
Reference  1.  Through  use  of  selective  computer  runs  of  modified  REPSIL  code 
it  was  possible  to  evaluate  the  relative  importance  of  each  of  these  effects 
on  the  overall  response. 

The  3003  aluminum  alloy  used  in  fabricating  the  beam  specimens  was  deter¬ 
mined  by  Hopkinson  bar  experiments  to  be  significantly  more  strain-rate  sen¬ 
sitive  than  most  structural  aluminum  alloys.  The  use  of  the  hyperbolic  tan¬ 
gent  function  strain-rate  model  introduced  in  this  report  is  believed  prefer¬ 
able  to  the  widely  used  power  law  model  due  to  its  more  realistic  behavior 
for  small  strain  rates.  Inclusion  of  strain-rate  effects  in  the  response 
calculations  for  Shot  5  produced  only  a  small  change  in  beam  deflection  due 
to  the  rather  small  strain-rates  which  were  associated  with  this  shot.  One 
should  not  conclude  that  strain-rate  effects  are  always  insignificant;  for 
problems  where  strain-rates  in  the  saturation  region  occur  for  an  appreciable 
portion  of  the  response,  resisting  stress  increases  of  40%  and  7R%  would  he 
available  for  Fits  I  and  II,  respectively. 

Although  the  intended  simply  supported  end  conditions  for  the  beam  tests 
were  not  achieved  experimentally  the  boundaries  were  certainly  not  completely 
fixed.  The  augmentation  of  REPSIL  boundary  conditions  to  account  for  support 
friction  and  rotatory  inertia  permitted  the  determination  that  the  observed 
response  was  primarily  due  to  large  frictional  forces  at  the  boundaries.  The 
relatively  small  effect  of  the  rotatory  inertia  of  the  rocker  shafts  did  alter 


the  mathematical  character  of  the  problem  and  provided  a  smoothing  effect  on 
the  response  calculations. 


A  motivating  consideration  for  performing  the  cited  beam  tests  had  been 
to  establish  the  existence  of  large  transient  shears  near  the  ends  of  blast- 
loaded  simply  supported  beeuns.  The  presence  of  such  shears  could  necessitate 
a  revision  of  statics-based  design  procedures  for  beams  and  slabs.  Calcula¬ 
tions  for  all  boundary  conditions  used  in  this  investigation  indicate  that 
such  shears  are  present  and  this  is  consistent  with  the  test  data.  However, 
unlike  the  simply  supported  beam  problem  in  which  bending  is  insignificant 
where  the  shears  are  large,  for  fixed  or  "locking"  ends  the  region  near  the 
ends  of  the  beam  is  subjected  simultaneously  to  large  bending  and  shearing 
stresses.  Therefore,  design  procedures  for  beams  with  constrained  ends  must 
consider  the  dynamic  combined  stress  states  which  will  be  present. 


''w.  A.  Kpryinn,  "Shr’ar  Stress  in  One-Way  Slabs  Sub{jected  to  Blast  Load,"  Civil 
Enyineerirui  Laboratory  Tochnioal  Memorandum  No.  M-51-77-10 ,  August  1977. 
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